Gravitational- Wave by Binary OJ287 in 3.5 PN Approximation 
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We compute the gravitational radiation ol the supermassive binary black hole OJ287. By fitting 
the data of its recent seven outbursts, we obtain the orbital motion up to 3.5PN (Post-Newtonian) 
order, the energy and angular momentum fluxes and the waveform up to the 2PN. It is found that 
the 1PN term of the energy flux, has an opposite sign to the Newtonian flux, and the 2PN term has 
an opposite sign to the 1PN one. The same pattern is also found for the angular momentum fluxes. 
The total flux is reduced by 30.8% on average from the Newtonian flux for OJ287, resulting in a 
smaller orbital decay rate than that of 2.5PN calculation. Consisting with this, it is checked that, in 
the sequence of non-dissipative PN forces, each term of i th order (dissipative or non-dissipative) has 
a sign opposite to the (i — l) th order. The origin of this characteristic is traced to the appearance 
of the non-diagonal metric component ga in PN approximation. This feature will have profound 
impact on estimation of gravitational waves from binary systems. 



Introduction Blasar OJ287 at a redshift z — 0.306 has 
long been observed continuously, and has been identified 
as a binary supermassive black hole system [lH5j with 
masses mi = 1.84 x 1O 1O M and m 2 = 1.46 x 1O 8 M . 
Having a very high orbital velocity v/c ~ (0.06 ~ 0.26), it 
serves as a test platform for general relativity up to 3PN 
Q, mainly via the precession effect on the orbit. Nev- 
ertheless, OJ287 is also supermassive, its gravitational 
radiation is very strong with amplitude h ~ 10~ 15 as 
we shall see, and is a proper target of the low-frequency 
gravitational wave detectors. In this paper we compute 
its orbital motion up to 3.5PN, i.e, including the second 
order of dissipative term. Then we obtain for the first 
time its gravitational waveform up to 2PN, and demon- 
strate the energy and angular momentum fluxes up to 
2PN correction to the Newtonian fluxes. The result is 
comprehensive and can be used for waveform templates 
for detectors. 




FIG. 1: After ~200 periods, 3.5PN orbit (solid) surpasses 
that of 2.5PN (dash), both starting with the same initial po- 
sition. 
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Orbital Motion to 3.5PN Since mi 3> mj, the con- 
figuration of the binary is modelled as such that the sec- 
ondary is orbiting around the primary, the accretion disk 
of the primary is perpendicular to the orbital plane . 
We do not include the possible spin effects for simplicity 
Q- The 3.5PN equation of motion in the center-of-mass 
frame is given by [8l4l4|. 



dv 
~dt 



Gm x 

[(1 + ^4iPN + ^2PN + ^2.5PN + ^3PN + ^3.5PN)e r 
(£>1PN + -B2PN + £?2.5PN + -B3PN + -B 3 . 5 pn)v], (1) 



where m = m% + m 2 , v is the relative velocity, e r is the 
radial unit vector, Aip^ and Bjpn are the PN correc- 
tions of the i order. In the polar coordinate (r, $) on the 
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orbital plane, the solution is obtained straightforwardly 
15], given an initial condition: ro, ro, $0, $o- By fitting 
with the observational data of the last seven outbursts of 
OJ287 B-fl, and taking into account of the time delay 
effect [15j , using the method of least squares for the differ- 
ence between the predicted and observed times of seven 
outbursts, we arrived at the optimal orbital parameters 



mi 
m 2 

period P 

ellipticity e 

semi — major axis a 

precession rate \& 



1.84 x 10 lo 7W Q 

I. 46 x 1O 8 M 

II. 9248 years 
0.6693 
0.0535 pc 
238.35° at 1971.48 
32.6359°/period 



The overall behavior of orbital motion is that the sec- 
ondary is orbiting around with a very large ^ due to 
the PN forces, the orbit is continuously shrinking and 
becomes more circular in the course of time, due to the 
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FIG. 2: The maximum mutual separation r max (left) and 
speed Umax (right) in 3.5PN are decreasing more slowly than 
those in 2.5PN, respectively. 



FIG. 3: Left: A2.5PN and A3.5PN have opposite signs at any 
instance. Right: -B2.5PN > and B3.5PN < 0. 



reaction of GW radiation. We have checked that the 
3.5PN orbital motion fits the data better than the 2.5PN 
one by the least squares 0, @. The detail has been 
given in Ref.[l5j]. As Fig Q] shows, the 3.5PN orbit has 
a slightly greater precession and gradually surpasses the 
2.5PN orbit. More interestingly, as shown in Fig(2J the 
mutual separation and velocity in the 3.5PN calculation 
are shrinking more slowly than those in the 2.5PN one, 
indicating that the binary described by 2.5PN must ra- 
diate away more gravitational energy than the 3.5PN de- 
scription. Now we investigate this feature in detail. 
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The Signs of PN Forces The first dissipative terms 
are at 2.5PN ones d[l3|, 
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with v = m\m2jm 2 . They cause damping of the or- 
bit and give rise to the Newtonian quadrupole radiation. 
The next order dissipative terms are 3.5PN, 



Without these four dissipative terms, as has been 
checked, the 3PN orbit will precess only and the radius 
will not shrink, so that the energy and angular momen- 
tum are conserved. 

FigEl demonstrates that, A2.5PN and A3.5PN always 
have opposite signs at any instance of time, and the am- 
plitude of A3.5PN is ~ 0.3 that of A 2 .5pn- So is with 
-B2.5PN and -B3.5PN- Therefore, the 3.5PN force is just in 
opposite direction to the 2.5PN one, the 3.5PN descrip- 
tion yields less orbital damping and less gravitational ra- 
diation than the 2.5PN one. Similarly, we find that in the 
sequence of the non-dissipative PN terms, Axpn, ^4 2 pn, 
^3PN, each term has a sign just opposite to its precedent 
term at any instance of time. The same alternating be- 
havior also occurs for the sequence of -Bipn, -B2PN, -Bspn- 
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The Energy and Angular Momentum Fluxes Since 
the dissipative 3.5PN forces are opposite to the 2.5PN, 
the work done (A3.5PNe r + -B3.5PNV) • dr is consequently 
opposite to (A 2 .5PNe r + -B 2 .5pnv) • dr. Thus one expects 
that, in the total energy flux [l6|-{l8[ 



dE_ 
~dt 



-En + -Eipn + -E2PN, 
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FIG. 4: The energy fluxes (in unit c 5 /G) E N > 0, E 1PN < 0, 

£<2PN > 0. 



the dominant Newtonian flux 
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induced by 2.5PN force has an opposite sign to 
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induced by 3.5PN force. In turn, for the same rea- 
son, £? 2 PN (see Refs. [TJ, [T3|) induced by 4.5PN force 
]19j has an opposite sign to Eipn 0- The energy 
fluxes of OJ287 are explicitly demonstrated in Fig|4]with 
i?N > 0, Eipn < 0, i?2PN > 0, which is in a com- 
plete agreement with the relative signs of 3.5PN and 
2.5PN forces. Similarly, the angular momentum fluxes 
H3 H El dJ/dt = j N + Jip N + j 2 PN of OJ287 is com- 
puted with a result that Jn > 0, Jipn < 0, and J 2 pn > 0. 
During each orbiting period the emission occurs mostly 
around the periastron where r is minimum and v attains 
the maximum. Integrated over a period, the negative 
Eipn reduces dE/dt by ~ 30%, and Jipn reduces dJ/dt 
by ~ 17%. This huge amount of reduction of radiated 
energy and angular momentum is due to the high or- 
bital speed of OJ287. Fig. [5] shows that the period P 
decreases more slowly in 3.5PN than in 2.5PN, a result 



that is explained by Eipn < and Jipn < 0. This pat- 
tern is consistent with the formula (4.29) of P and the 
formula_(4.31) of P in Ref.[JJ|, but different from (4.26) 
in Ref.[18|]. The lesson learned is that, 2.5PN calculations 
of the Newtonian fluxes En , Jn , and the heuristic rate P 
[2ll ] are over-estimations. A similar pattern has also been 
noticed in gravitational recoil due to linear momentum 
loss, where the Newtonia n, 1PN, and 2PN contributions 
change sign alternatingly |22H24j . 

The Waveform Up to 2PN, the traceless-transverse 
(TT) gravitational waveform is given by 16] 



R 



P 1/2 Q ij +PQ ij +P 3/2 Q 13 +P 2 Q ij ] TT 



where R is the distance from the binary to the observer, 
Q 13 is the dominant, quadrupole momentum part, and 
the remaining terms with superscripts denote the the ef- 
fective PN order (their explicit expressions are given by 
(6.11) in Ref . ) . Using the resulting orbital motion, we 
compute hj, T up to 2PN order. Since OJ287 is at cosmic 
distance with z ~ 0.306, the effect of cosmic expansion 
should be taken into account. The distance R is taken to 
be the luminosity distance: 



Dl{z) = + z) f dz ' == 



(9) 



in a spatially flat ACDM Universe with £7a = 0.685, 
H a = 67.3(km/s)/Mpc and fl m = 1 - Q A [25|, where H 
is the Hubble constant. For OJ287, R = D L (0.306) ~ 
1646.89 Mpc. The prefactor amplitude of hij is ~ 
8.42 x 10~ 15 . Let the observer's direction be given by 
n = sin #e x + cos 8e z , where e z is the normal vector and 
e x is along the semimajor axis on the orbital plane. Take 
an orthonormal basis (eg, e^, n). The two independent 
components of GW measured by the observer are [26} 



:;:' r = I( CO s 2 9h x 



h yy ), 



h"/ T = cos 9h xy . 



(10) 
(11) 



Fig|5] shows the wave form h^ T and hrf T using the 
3.5PN orbital result. The peak amplitude reaches \h 



13 



1 x 10~ 15 with a frequency / ~ 2.66 x 10~ 9 Hz, higher 
than the amplitude of relic gravitational waves in this fre- 
quency range [27J ■ Both can be the targets of pulsar tim- 
ing arrays, such as PPTA J2|, EPTA [29|, NANOGrav 
Ifl], FAST [HI and SKA [HJ But its frequency may be 
too low for LIGO and VIRGO 

The Origin of the Opposite Sign Now let us see the 

origin of the opposite sign occurring at each next or- 
der. Inspect that, aside the common factor, ^42.5PN in 
Eq.© contains the negative factor - if ^?), 

while the dominant terms of A3.5PN in Eq.(j3]), such as 



f 4 , v 2 ^, 



■2Gm j Gm\2 



(^T 2 ) , are positive. Thus A3.5PN 
has an opposite sign to A2.5pn- Similarly, inspection of 
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FIG. 5: The period P in 3.5PN decreases at a slower rate 
than in 2.5PN. 



0.10- 






























0.05- 
































































O.OO- 
































































-0.05- 
































































-0.10- 






























































-0.15- 

































0.15-, 


































0.10- 




































































0.05 - 




































































O.OO - 






















































s 


















-0.05 - 








































































0.10- 


























































5.5PN 








0.1 5 - 



































t/(Gm,/c 3 ) 



FIG. 6: The waveform h e TT (upper) and hrf T (lower) in 
3.5PN. The observer is facing the orbital plane with the angle 
= 0. 



Eq.© tells that -B2.5PN > 0, while the dominant terms of 



-B3.5PN in Eq.©, such as v 4 , r 
all have large negative coefficients, so -B3.5PN < 0. 

This feature can be inferred generally from the PN 
approximation, say, in the method of the PN iterative 
metric parameterized by retarded potentials 14]. The 



acceleration is determined by a 1 - 
Eq.(2.13) in Ref.0]), where 

F 1 ~ diV + ^(-VdiV + ~d t V v 2 



di 



(P l - v l ) (see 



-AB i V s v i ) + 0{l/c% 
AV l ) + 0(l/c 4 ), 



with V ~ Gm/r ~ v 2 and Vj ~ v 3 V to leading or- 
der. For each order, the dominant contribution from diV 
is followed by the next order contribution \{—Vd%V + 
%diVv 2 - AdiVjvi) ~ -l^v 2 diV, and -£(P l - v*) ~ 
-iUb 2 v l - Vv<) - pib 2 ^ - M^iV. The 
sum of these two terms is ~ — -^2v 2 diV having an op- 
posite sign to diV. Thus, for each PN order of a 1 , the 
next order has an opposite sign. The cause is traced to 
idiVj t)i) in F l , which originates from the metric 
component g oi — —Wij<? + 0(l/c 5 ), i.e., the shift [ilj] . 

Conclusion We have studied gravitational radiation of 
the binary black hole OJ287, apply the Post-Newtonian 
approximation up to 3.5PN order in the orbital motion. 
By explicit computations, we have found that, at any 
instance of time, the energy flux E\pn induced by the 
3.5 force is negative, opposite to the Newtonian flux En 
induced by the 2.5PN force, and, consistent to this, the 
dissipative 3.5PN force is always opposite to the 2.5PN 
one. For OJ287 this reduction of energy flux is as high 
as ~ 30%. Therefore, for binaries with high orbiting 
speed, this effect is significant and must be taken into 
account in evaluation of radiation of GW. We have also 
demonstrated that, for non-dissipative PN forces, each 
PN order has a sign opposite to the precedent PN or- 
der. The origin this characteristic has been traced to 
the non-diagonal metric component goi in PN approxi- 
mation. With a peak amplitude \hij\ ~ 10~ 15 and low 
frequencies / - 2.7 x 1(T 9 Hz, the GW from OJ287 is 
an ideal target for pulsar timing array detectors. 
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